Extended Pairing Model for Heavy Nuclei 
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We study binding energies in three isotopic chains ( 



3 Sn, 152 " 181 Yb, and 



J Pb) using 



the extended pairing model with Nilsson single-particle energies. The exactly solvable nature of the 
model means that the pairing strength G(A) required to reproduce the experimental binding energies 
can be determined uniquely. The valence space consists of the neutron single-particle levels between 
two closed shells corresponding to the magic numbers 50-82 and 82-126. In all three isotopic chains, 
\og(G(A)) has a smooth quadratic behavior for even as well as odd nucleon numbers A; \og(G(A)) 
for even and odd A are very similar. Remarkably, G(A) for all the Pb isotopes can be also described 
by a two parameter expression that is inversely proportional to the dimensionality of the model 
space. 

PACS numbers: 21.10.Dr,21.60.Cs, 03.65.Fd,71.10.Li,74.20.Rp,27.70.+q,02.60.Cb 



In a previous paper we introduced the extended pair- 
ing model, explored its properties, and applied it to some 
well-deformed nuclei [lj. Metallic clusters of nano-scale 
size may provide another physical system where one can 
test and study the applicability of the extended pair- 
ing interaction 0. In this paper we exploit the exactly 
solvable nature of the model to explore its applicability 
in reproducing experimental nuclear binding energy of 
three distinct isotopic chains: 100 ~ 130 Sn, 152_181 Yb, and 
181 ~ 202 Pb. The results suggest that the model is useful 
for its simplicity, its exact solvability, and its ability to 
track and predict experimental binding energies of long 
sequences of nuclei. 

For each isotopic chain we take the binding energy of 
a closed neutron shell to be our zero-energy reference 
point. Each such closed neutron shell nucleus ( 100 Sn, 
152 Yb, and 208 Pb) and its odd- A neighbor ( 101 Sn, 153 Yb, 
and 207 Pb) are assumed to be well described by the in- 
dependent particle model with Nilsson single-particle en- 
ergies; thus no extended pairing interaction terms are 
needed for these nuclei. The energy scale of the Nilsson 
single-particle energies is set so that the binding ener- 
gies of 101 Sn, 153 Yb, and 207 Pb are reproduced by the 
independent particle model. For all the other nuclei we 
solve for the pairing strength G(A) that reproduces the 
experimental binding energies exactly within the chosen 
model space. The valence model space consists of the 
neutron single-particle levels between two closed shells 
corresponding to the magic numbers 50-82 and 82-126. 
The structure of the model space is reflected in the values 
of G(A). In particular, in all the cases studied log(G(A)) 
has a smooth quadratic behavior for even and odd A 
with a minimum in the middle of the model space where 
the size of the space is a maximal; log(G(A)) for even A 



and odd A are very similar which suggests that more de- 
tailed shell- model analyses may result in the same G(A) 
functional form for even A and odd A isotopes. In par- 
ticular, for the Pb isotopes the even and odd A nuclei 
are described by a single two-parameter expression for 
the pairing strength G(A) that is inversely proportional 
to the dimensionality of the valence model space. 

To set the stage and establish the notation we begin 
with a brief review of the underpinnings of the extended 
pairing model. The standard Nilsson plus pairing Hamil- 
tonian is given by 
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where p is the total number of single-particle levels con- 
sidered, tj are single-particle energies, G is the overall 
pairing strength (G > ), rij — c ]| c j'T+ c ]i c j1 i s ^ ne num_ 
ber operator for the j-th single-particle level, af = ct-.ct, 
{fli = (a^Y = Cj|Cii-) are pair creation (annihilation) 
operators where cl^ (cji) creates a fermion in the j-th 
single-particle level. The up and down arrows refer to 
time-reversed states. Since each Nilsson level can only 
be occupied by one pair due to the Pauli Exclusion Prin- 
ciple, the operators af, a«, and rii form a hard-core boson 
algebra: 



[a i ,af]=8 ij {l-n i ), [a+ , af] = = (a+f . (2) 

As a generalization of Q}, we consider the following 
extended pairing Hamiltonian [l|: 
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Besides the usual single-particle terms and the standard 
pairing interaction this interaction includes many- 
pair interaction terms which connect configurations that 
differ by more than one pair. 

The main advantage of the extended pairing model 
is that it is exactly solvable pj. It is easy to see that 
any term of the form af ■ ■ • at that forms eigenstates of 
© should enter with different indices i ^ • • • ^ j. Let 
' • • ,jm) denotes a pairing vacuum state that satis- 
fies Oi|ji, • • • ,j m ) = for 1 < i < p and i ^ j a , where 
iii 3i i • • ■ 7 jm indicate those m levels that are occupied 
by single nucleons. Any state that is occupied by a sin- 
gle nucleon is blocked to the hard-core bosons due to the 
Pauli exclusion principle. This means that the space of 
all possible configurations decomposes in a direct sum of 
orthogonal sub-spaces that are invariant under the ac- 
tion of the Hamiltonian and are labeled by the positions 
of the unpaired nucleons. Thus a fc-pair eigenstates of 
(0 has the form 



allowed. It is always assumed that the level indices 
ji j J2 1 • • • j jm should be excluded from the summation 
in Q. Since the general formalism is similar, we will 
focus on the seniority zero case (no unpaired particles). 
In a manner that is similar to the results given by the 

Bethe ansatz, the expansion coefficient Cj^.., in Q 
can be written as 0: 
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where y^' is a number that is to be determined. To solve 
(0 using and (J5J one directly applies Hamiltonian 
© on the fc-pair state (QJ. Using the hard-core boson 
algebraic relations given by (|2J), one can determine the 
action of the mean-field part of the Hamiltonian J3J) : 



E' 



(6) 



|fc;C;ii, • • • J m ) = c h--ik a h ' ' ' a t I-? 1 ' ' ' ' >3m), 

l<ii<--<ik<p 

(4) 

where C^J ^ are expansion coefficients that need to 
be determined, and the strict ordering to the indices 
*ij *2> * " " >ik reminds us that double occupation is not and for the extended pairing part of the Hamiltonian J3J): 
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By combining © and J7J), the /c-pair excitation energies 
of @ are given by 

EP = -^-G(k-l), (8) 



and the undetermined variable is given by 
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The additional quantum number £ can now be under- 
stood as the £-th solution of 10. Similar results for 



many broken-pair systems can be derived by using this 
approach except that the indexes j s of the level occu- 
pied by the single nucleons should be excluded from the 
summation in Q and the single-particle energy term tj s 
from the first part of should be added to the total 
eigenenergy. 

By comparing JSJl and @ to the exact solutions of 

the Heisenberg algebraic Hamiltonian with a one-body 

interaction [4j, one can regard the operator product 

af af ■ ■ ■ af in H as a 'grand' boson. The Corre- 
al 22 * i » G 

sponding 'single-particle energy' of the 'grand' boson 
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The eigenstates (@J are not 
finalized once the coefficients 
are known. The eigenstates |(2J with different 



is Ei 1 i 2 ,_i k — 

normalized, but can be normalized once the coefficients 



roots given by © are mutually orthogonal since they 
correspond to eigenstates with different eigenvalues. 

The k coupled non-linear equations of the standard 
pairing model [f| are difficult to solve numerically, espe- 
cially when the number of pairs k and number of levels p 

dis- 



are large. Specifically, there should be K 



(p-fc)!fc! 



tinct roots, which can be a very large number for an entire 
deformed major shell. While there have been major ad- 
vances in methods for solving the associated Richardson 
equations, for example [6(, the theory is limited because 
of the coupled non- linear nature of the equations. In con- 
trast to this, for the extended pairing model there is but 
one variable and one equation @; so a relatively 
simple Mathematica code, for example, can be used to 
solve for the roots. 

Our Nilsson plus extended pairing model uses single- 
particle energies of each nucleus as calculated within the 
Nilsson deformed shell model with experimentally evalu- 
ated deformation parameters Q. Experimental binding 
energies are taken from reference Q. The theoretical 
binding energies are calculated relative to a particular 
core. We use 152 Yb, 100 Sn, and 208 Pb as cores in our 
calculations. While there are changes in the binding en- 
ergy of the core since the corresponding Nilsson levels 
change as a function of the deformation, our results in- 
dicate that such core affects are an order of magnitude 
smaller than the overall even-odd staggering in the bind- 
ing energy. From the binding energy of a nucleus next to 
the core, we calculate an overall energy scale for the Nils- 
son single-particle energies. For an even number of neu- 
trons, we consider only pairs of particles (hard bosons). 
For an odd number of neutrons, we apply Pauli block- 
ing of the Fermi level of the last unpaired fermion and 
consider the remaining fermions as if they were an even 
A fermion system. By using JSJl and ©, values of G 
are calculated so that the experimental and theoretical 
binding energy match exactly. With the help of JSJ and 
@ and some algebra, one can see that for a given single 
particle energies there is an upper limit to the value of 
the binding energy for which a physically meaningful ex- 
act solution can be constructed. This upper value of the 
binding energy for each nucleus is given by the energy of 
the lowest "grand boson", namely, E g b = Ylu=i 

We now turn to our results for the binding energies and 
log(G) values within the Nilsson plus extended pairing 
model for three isotopic chains: 100_130 Sn, 152 ~ 181 Yb, 
and 181 ~ 202 Pb. Figure 2] shows the results for the 
ioo-i30g n i go t p es Calculations for the pairing strength 
G were carried out for the 102_130 Sn isotopes within the 
50-82 neutron shell. Once this was done, a quadratic 
polynomial fit of the log(G) values for even and odd A 
was determined. By doing this we were able to fit two 
sets of 14 data points with two 3 parameter expressions. 
Overall the results are very good for the lower and middle 
part of the 50-82 neutron shell. Even though a particle- 




FIG. 1: The solid line gives the binding energies of the 
Sn isotopes relative to that of the 100 Sn core. The single- 
particle energy scale is set by the binding energy of 101 Sn. 
The inset shows the fit to values of G that reproduces the 
experimental data exactly. The two fitting functions are: 
log(G(yl)) = 365.0584 - 6.4836A + 0.0284.4 2 for even values 
of A and log(G(A)) = 398.2277- 7.0349A + 0.0307A 2 for odd 
values of A. The Nilsson BE energy is the lowest energy of 
the non-interacting system. 



hole symmetry can be seen from the log(G) inset, there 
is a discrepancy between the two as one moves towards 
the upper part of the shell. This discrepancy is due to 
the shell closure, the single-particle level structure, and 
the Pauli blocking for odd A nuclei. 
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FIG. 2: The solid line gives the binding energies of Yb isotopes 
relative to that of the 152 Yb core. The single-particle energy 
scale is set from the binding energy of 153 Yb. The inset shows 
the fit to values of G that reproduce exactly the experimental 
data. The two fitting functions are: log(G(yl)) = 662.2247 - 
7.7912,4 + 0.0226.4 2 for even values of A and log(G(yl)) = 
716.3279 - 8.4049,4 + 0.0244.4 2 for odd values of A. The 
Nilsson BE energy is the lowest energy of the non-interacting 
system. 

Figure |3 shows our results for the 154 - 161 Yb iso- 
topes. We test the predictive power of the model on 
the 172 - 177 Yb isotopes. In this case calculations of 
the pairing strength G are carried out only for the 
154-171, 178-181 yb isotopes and not for the ^-myb iso _ 
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topes that are in the middle of the model space and thus 
are more computationally involved. This is a fit of two 
sets of 11 data points with two 3 parameter expressions. 
Then from the obtained quadratic polynomial fit to the 
log(G) values we calculate the theoretical values of the 
binding energy for these nuclei as shown in Figure El 
This prediction is very good when compared to the ex- 
periment. Thus, based on experimental data of the nuclei 
in the upper and lower part of the shell and a log(G) fit to 
this data we can make a reasonable estimate for the mid- 
shell nuclei. Therefore, the model has a good predictive 
power. 




FIG. 3: The solid line gives the binding energy for the Pb iso- 
topes relative to the 208 Pb nucleus. The insets show the fit to 
the values of G that reproduce exactly the experimental data 
using a 164 Pb core. The lower inset shows the two fitting func- 
tions: \og(G(A)) = 382.3502 - 4.1375^1 + O.Olllyl 2 for even 
values of A and log(G(A)) = 391.6113 - 4.2374A + 0.01144 2 
for odd values of A. The upper inset shows a fit to G(A) 
that is inversely proportional to the size of the model space, 
(dim( J 4)), that is valid for even as well as odd values of A: 



G(A) = 366.7702/ dim(A) 



The Nilsson BE energy is 



the lowest energy of the non-interacting system. 

The next Figure shows results for the 181 - 202 p D iso- 
topes that were studied in the same way as the Sn and 
Yb isotopes. In the calculations for these Pb isotopes, 
however, the binding energy that was used is relative to 
that for 208 Pb which was set to zero, but the core nu- 
cleus was chosen to be 164 Pb. Note that for the Yb and 
Sn isotopes the core nucleus was also the binding en- 
ergy reference nucleus ( 100 Sn and 152 Yb). In contrast, 
the calculation for the Pb-isotopes was different because 
the core nucleus ( 164 Pb) and the binding energy refer- 
ence nucleus ( 208 Pb) are different. We again have good 
quadratic fit to log(G) as function of A. 

The fact that there is a correlation between the pairing 
strength G and the size of the model space reflected in the 
minimum of G that is at the maximal model space dimen- 
sion prompted us to study G(A) as function of the model 
space dimension dim(A). In this respect, a remarkable 
results is shown in Figure |3J In this case the pairing 



strength G(A) for all the 21 nuclei (A=181 - 202) was fit 
by a two parameter function G(A) — a/[dim(A)]P with 
the values of the parameters taken to be a — 366.7702 
and (3 = 0.9972. This function is inversely proportional 
to the dimensionality of the model space dim(^4). Unfor- 
tunately, this is not the case for the other two isotopic 
chains. For example, in Figure 01 a log — log plot is shown 
for the Sn nuclei. 
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FIG. 4: Results for log(G) versa log(dim) for the Sn nuclei. 

It is unlikely that the linear relation between log(G) 
and log(dim(j4)) for the Pb isotopes is due to the dif- 
ference between the core and binding energy reference 
nucleus. What is more likely is that this is due to the 
fact that in these cases the model space dimensions are 
sufficiently large to result in a limiting form for the effec- 
tive pairing strength. Another reason might be that the 
single-particle energies and the Pauli blocking mechanism 
are such that the even and odd parabolas of log(G(A)) 
produce a linear log — log structure. 

In the light of the above, it seems that the next step 
should be a study that tracks the results as a function 
of the increasing size of the model space to confirm or 
refute the log — log relation. Such a study could also ad- 
dress other questions such as the effect of the core bind- 
ing energy as a function of the deformation that is used 
in the Nilsson model to derive the single-particle ener- 
gies. Using a Woods-Saxon potential or other methods 
to generate more realistic single-particle energies is an- 
other opportunity for further studies. 

In conclusion, we have studied binding energies of nu- 
clei in three isotopic chains: 100_130 Sn, 152 - 181 Yb, and 
181 ~ 202 Pb within the recently proposed extended pairing 
model by using Nilsson single-particle energies as the 
input mean-field energies. Overall, the results suggest 
that the model is applicable to well-deformed nuclei if 
the pairing strength is allowed to change as a (smooth) 
function of the nucleon number A. The remarkably simi- 
lar behavior of log(G) for even and odd A values seems to 
suggest that there may be a single log(G) function that 
bifurcates into an even- A and an odd- A branch when the 
fermion dynamics is restricted to hard boson pairs only 
and Pauli blocking is applied to exclude levels populated 
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by the unpaired fermion. 
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